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Abstract. Our fundamental lack of understanding of the acceleration of the 
Universe suggests that we consider a kinematic description. The simplest for- 
malism involves the third derivative of the scale factor through a jerk parameter. 
A new approach is presented for describing the results of astronomical observa- 
tions in terms of the contemporary jerk parameter and this is related to the 
equation of state approach. Simple perturbative expansions about ACDM are 
given. 



1. Approaches to Dark Energy 

The Universe is spatially flat, accelerating and has a subcritical matter density. 
The evidence for this comes from: 

• observation of acoustic peaks in the microwave background radiation spec- 
trum which effectively measures the angular diameter distance dA to the 
last scattering surface (e.g. Spergel et al. 2003). 

• magnitudes of Type la supernovae which provide the luminosity distance 
di as a function of redshift (e.g. Riess et al. 2001). 

• X-ray observations of rich clusters of galaxies which measure the mean 
density of dark matter and also, with less confidence, the distance-redshift 
relation (e.g. Allen et al. 2004). 

Although the possibility that the Universe contains a repulsive entity coun- 
teracting the attractive effects of gravity was entertained by Einstein in the 
first, relativistic cosmology, the discovery of acceleration was mostly a surprise. 
There have been several approaches to rationalizing this discovery from a con- 
temporary perspective. The simplest, and probably still the favorite description, 
is that the field equation is supplemented by a term proportional to the met- 
ric tensor, which guarantees covariance. The coefficient of proportionality (the 
cosmological constant) is the single free parameter in the theory. 

If we regard this term as an augmentation of the stress energy tensor of cold 
matter and adopt a fluid description, then the enthalpy must vanish and the pres- 
sure is the negative of the energy density (roughly 0.7 nJ m~^). Alternatively, 
we can regard this term as a necessary component of the geometrical response 
of a spacetime manifold to a material source that only becomes detectable on 
cosmological scales. Either choice leads to the flat "ACDM" cosmology in which 
the cosmic time varies with the scale factor a = 1/(1-1-2;) according to the 
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solution of Lemaitre (1927) equation 
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(Bondi 1952), where Q^n is the density parameter of matter (p = Qpc where 
the critical density pc = SHq/SttG) and a is normalized to unity at the present 
time-^. Following Spcrgcl et al. (2003), we adopt Hq = 72 km s^^ Mpc"-*-. 

In order to test this description, it is necessary to introduce a framework 
which embraces a larger class of models labeled by measurable parameters. The 
overwhelmingly favourite choice of model for "dark energy" is to suppose that 
the acceleration is caused by a fluid which expands according to a modified 
adiabatic law, so that PV^~^^ is equal to a (negative) constant (Steinhardt 1997; 
Turner & White 1997) . (This is reminiscent of a polytropic description of a star 
with the important distinction that it is the same material that undergoes the 
change as the Universe expands as opposed to matter with different thermal 
histories as is the case with a star.) The parameter w (= —1 under ACDM) can 
be allowed to vary but wc shall treat it as constant for the moment. Given this 
prescription a satisfies a Friedmann differential equation (neglecting curvature) 



This approach can be confusing because, if we were to treat w as a pure 
number then small scale perturbations to the dark energy fluid would grow ex- 
ponentially with time for w < 0. This is clearly untenable. Instead w is usually 
associated with a scalar field theory with a more complex rule for handling 
perturbations. Under these circumstances, it would seem preferable to parame- 
terize the field theory itself which will lead to a different functional departure of 
the models from ACDM. An analogy is that one invariably solves problems in 
electromagnetic theory working from Maxwell's field equations rather than his 
stress tensor. In practice small scale dark energy perturbations propagate with 
the speed of light in most, though not all, models (e.g. Liu 2004). 

A huge variety of scalar field models have been discussed (e.g. Caldwell et.al. 
1998; Peebles & Ratra 1988). However, there is no guarantee that the accelera- 
tion of the Universe has anything to do with a scalar field. Other explanations in- 
clude the notion that the Universe has extra dimensions (e.g. Dvali Sz Tye 1999). 
The scalar field is probably the best bet although it should be remembered that 
no such fundamental fields have been detected experimentally as yet. Electro- 
magnetic and gravitational fields were thought to be scalar, but turned out not 
to be. To make this point, consider a blind astrophysicist practising her craft 
during the radiation era. She might have measured that the average spacing of 
protons increased as the square root of time and have inferred that the expan- 
sion was driven by a scalar field with an exponential potential. She would have 
been hopelessly wrong. 



2H^V{a) = Hi ^rna-^ + (1 - n^)a-(i+3"') 
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^Eddington (1935) strongly advocated a ACDM model on the grounds that it would prolong the 
life of the Universe and that the existence of a cosniological constant was necessary to explain 
the disparity between the atomic and cosmological scales. 
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2. Kinematic Approach 
2.1. Jerk Parameter 

In view of the huge diversity of dynamical models on offer, there may be some 
merit in the observational cosmologist taking an unprejudiced kincmatical ap- 
proach. This has been historically productive. Galileo, who originally argued 
theoretically that the speed of a falling body increased in proportion to distance, 
determined the true law by empirical, kinematic measurement (Dugas 1955). 
Modern cosmology began with the discovery of the expansion of the universe, 
quantified by the Hubble parameter H{t) = d/a. In more recent times the 
pursuit of observational cosmology was dominated by the measurement of the 
contemporary deceleration parameter go where q{t) = —da/a^ (Sandage 1975). 
In terms of the scale factor and Hubble parameter q = —d In aH/d In a. 

As the Universe was once decelerating and is now accelerating, it is use- 
ful to consider the third derivative of a. A convenient quantity to use is the 
dimcnsionless jerk (Blandford 2004) 



with its current value denoted by jo- One reason why this is convenient is that 
j = 1 for ACDM, the baseline model about which we are perturbing. In terms 
of the Hubble parameter j = (a^ H^)" /2H'^ . Different formalisms involving 
the third derivative have been discussed independently by (Sahni et al. 2003; 
Visser 2004). 

2.2. Dynamical considerations 

Although, we eschew dynamics, j can be interpreted in fluid terms, through the 
easily derived relation j = I — AttP/H^ = 1 + {a^p')' /2pa'^ , where P is the pres- 
sure, p the total energy density while we ignore spatial curvature (analagously 
—q = l + ap'/2p). The expressions involving H rather than p above are valid 
even considering curvature. For a slow rolling scalar field (f) and potential U {(j)) 
we get i « 1 - 4TrU'{(l))^/3H^. 

2.3. Hubble Parameter 

We solve Eq. 3 by rewriting it in the form 



where V{a) = — d^/2i?o = —H^{a)a^ /2Hq and prime denotes differentiation 
with respect to a. This equation has power law solutions for constant j, but the 
only one that is of interest is the ACDM solution with j = 1. This is because we 
seek solutions that approximate the Einstein-De Sitter case for early times. We 
therefore take for our simplest, one parameter extension of ACDM, the variation 
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j(a) = l+/a = l + (jo - l)a 



(5) 



In general, if we specify an expression for j{a) and require that the solution 
match an Einstein-De Sitter solution at early time, then we can solve for the 
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Figure 1. The effective potential V{a)^ and the derivative din V{a) /df , dw, 
evaluated at / = (w = — 1). 



scale factor, a{t), the Hubble parameter, H(t), and the deceleration parameter, 
q{t) for all time. (Specifically all we need to know is j{a) and lim.a^o{aa'^ / Hq) .) 
As is the case with the w formalism, it is possible to introduce extra parameters 
through a Taylor expansion j(a) = 1 + j'a + j"a^/2 + . . . However, we shall not 
explore this generalization here, and instead restrict ourselves to a jerk that is 
linear in the scale factor. Eq. 4 has the solution 

V{a) = -I \cifi{8aj')/a + C2f2{8aj')a^] (6) 



where 



= x^/^K3{x^/^)/8 = l-x/8 + x'^/64 + ... 



(7) 

and ci^2 are constants. 
2.4. j vs. w 

In order to compare the j and w prescriptions, we must compare world models 
computed adopting similar boundary conditions. We choose to fix the co-moving 
distance Hod{ar) = Hq dt/a to the recombination surface(ar ~ 0). Note that 
the distance to the recombination surface, HQd{0) is in practice measured quite 
accurately by the microwave background fluctuation spectra. In everything that 
follows, we shall assume HQd{0) = 3.4. This condition allows us to solve for 
as a function of j' (or w). If we linearize in j' and assume that matter evolves 
in the usual way then in the j model, ci = flm and C2 can be chosen such 



Cosmokinetics 



5 




0.75 



Figure 2. The kinematic jerk parameter j(a) shown for three different val- 
ues of the dark energy equation of state parameter w. 
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-1/2. The linearized expressions are shown below for illustrative 
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(8) 



In Fig. 1 we plot the exact solution for V{a) along with the derivative din V{a) /dj' , dw 
evaluated at j' = {w = —1) to show the difference between the w and j for- 
malisms. 

In general dlnX/dj' characterizes the deviation of an observable X from 
the ACDM model. This can be seen by considering the following expansion: 



X{a,f)^Xia,0){l+j 



,dlnX 



dj' 



j'=Oj 



(9) 



Note that this provides a simple route to the generation of fitting formulae. 
Substituting w for j' and expanding about the ACDM value of it; = —1 gives 
the corresponding expression in the dynamical formalism. 
We can also compute j{a) for a given w model. 



J{a) 
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This relationship is plotted in Fig. 2. 
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2.5. Distance 

Many techniques for exploring the expansion of the Universe utiUze distance 
measurement through the angular diameter dA = ad or luminosity distance 
= a~^d or the comoving volume V = Aird^ /?,. For ACDM, 

where B[a,b;x] = Jq dtt"'~^{l — t)^~^ is the incomplete Beta function. The 
distances can be computed in the perturbed models with the the boundary 
condition HQd{0) = 3.4. We contrast the j and w formalisms in Fig. 3 where 
we plot Hod and dlnHod/dj' ,dw evaluated at / = 0{w = —1) as functions 
of a. The fractional changes in the co-moving volume are three times these 
functions. A simple fitting formula for H^d is given by combining Eq. 9 with 
dlnHod/df ^ En=iCna", where the coefficients c„ = {0.54,-1.1,0.81,-0.26}. 



2.6. Curvature 

The WMAP observations have demonstrated that the Universe is quite flat; 
any curvature can surely be treated as a perturbation and will be very hard 
to detect in the face of the above effects. We therefore contrast the devia- 
tion expected from pure curvature in the distance relation. Curvature with ra- 
dius R appears in standard Friedmann cosmology in two related places. The 
first is in the relation between the comoving radial coordinate r = J dt/a 
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Figure 4. Co-moving distance Hod(a), and the derivative din Hod{a)/dflk, 
evaluated at — 0,/ = 0{w = —1). 



and the distance d = r + Q^r^/G + . . ., and = where the sign in- 

dicates positive or negative spatial curvature (0^ negative is positive curva- 
ture). Note that this effect is purely geometric, making no reference to dy- 
namics. The second is in the dynamical equation, Eq. 2, which changes so 



that V 



+ (i-a 



m - Ofc)a-(^+3"') + VLk /2. The jerk parameter 

can be derived from j = 1 + {a^p')' 120? p by defining pk = 0,kPca~'^ and taking 
p ^ p + Pk- For the pure w model, jo = 1 — ilfc + 9(1 — Vim — VLk)w{l + w)/2. 
In Fig. 4 we show dlnHod/d^lk evaluated at fifc = 0,j' = 0{w = — 1). 



2.7. Time 

At present we do not have very accurate age measurements but it is still useful 
to compute the relative changes in the chronology induced by our models. We 
obtain the time under ACDM from Eq. 1. The times for the w,j modifications 
are obtained by integrating Eq. 8 subject to a third boundary condition that 
t{0) = 0. dlnHot/dj',dw for / = 0{w = -1) are plotted in Fig. 5. 



2.8. Growth of Perturbations 

The rate at which the perturbations that are observed directly in the microwave 
background fluctuation spectra grow under the action of Newtonian gravitational 
forces is quite sensitive to the expansion rate of the Universe. If we ignore 
radiation, pressure and the distinction between cold and baryonic matter, then 
a relative matter density perturbation on a flxed co-moving scale, denoted by 
(5, grows according to a second order differential equation which can be cast in 
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Figure 5. Cosmic time Hat{a), and the derivative dhiHot{a)/dj\dw, eval- 
uated at j' — Q{w — —1). 



Sturm-Liouville form 



(^.^vvy + 11^ = (12) 



We solve this differential equation evolving the linear growth transfer function 
for the w and j models from the epoch of recombination to the present day 
assuming 5{ar) = 1 and 5'{ar) = l/a^. The difference between the w and j 
models is shown using d\n6/dj' ,dw. 



3. Discussion 



The above calculations suffice to contrast the j and w formalisms. The procedure 
that is actually followed in fitting various data sets is to explore the likelihood of 
a set of cosmological parameters subject to prior assumptions about the nature 
of the solution. It is clear that sensitivity of a given analysis to a specific param- 
eter depends upon the source redshift distribution and such choices can either 
flatter or disparage a given approach! The two choices contrasted here are, in 
the absence of any compelling theoretical guidance, essentially arbitrary and it 
is not our purpose to advocate one over the other. In particular, it seems to be 
a mistake to suppose that because the w formulation produces maximal devia- 
tion from ACDM at a ~ 0.6 then observations should be concentrated around 
z ~ 0.7. Rather, we seek to raise awareness of alternative parameterizations 
of the cosmological model which may help increase understanding, and prevent 
misunderstanding, in the face of forthcoming data. 
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Of course, just as one should consider performing dark energy observations 
over the whole available redshift range, one should also consider a range of mod- 
els when analysing these data. Attempts to measure functions d{a) directly from 
the data via some necessary interpolation scheme form one such class of models. 
The much simpler jerk and dynamical parameterisations described above trade 
flexibility for ease of interpretation, but are as valuable and valid as analytical 
tools. The problem of which model provides the most appropriate description 
of these data has a well-defined solution that depends on both the goodness- 
of-fit (via the likelihood) and the model complexity (via the prior). The jerk 
parameterisation is a natural way to perturb the model that we know to fit the 
current data very well, and so may be expected to come with some sensible prior 
for the parameter j' . Indeed, noting that it is a straight line gradient, insisting 
on rotational invariance in the plotting plane leads to a Lorentzian prior pdf 
of unit width (Dose 2003). This distribution is broad enough to allow the data 
to speak for themselves, and featureless for as long as we ignore any dynamical 
considerations. The situation in the dynamical modelling scenario is somewhat 
different: there, the onus is on the theorists to provide prior distributions for 
physically-motivated parameters such as w, a task that has not been undertaken 
in any great detail as yet. 

Our main objective in this work has been to provide an alternative approach 
to the problem of characterising the cosmological world model in the presence 
of a poorly understood dark energy component. We may expect simple param- 
eterisations to be of great use in interpreting the data arriving from the many 
current and future dark energy experiments; this is especially true in the early 
stages when the signal-to-noise is low. The jerk formalism worked through here 
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has a number of advantages; we hope that it will be considered by observers 
as a useful parallel analysis tool, and as such we have provided an example fit- 
ting formula for the co-moving distance. A more comprehensive jerk toolkit, 
and a more detailed discussion of points made earlier in this talk, will be given 
elsewhere. 
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